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1. Introduction 



Recently, it was shown that the nonlinear evolution of anti-de Sitter (AdS) space might 
be unstable, and the energy of perturbations would be transformed to smaller and 
smaller scales like the turbulence energy cascades |], |j, and it is interesting to estab- 
lish the holographic turbulence using the AdS/CFT correspondence ||, f|. For example, 
in || |[ 0, H, H, the problems of turbulence have been studied via a holographic de- 
scription with gravity. Especially, using the fluid/gravity correspondence |1^, |TI| , it 
was suggested that the problem of Naiver-Stokes (NS) turbulence might be mapped 



to a problem in general relativity Jl2], [L3|], with different scales appearing in turbu- 
lent phenomena corresponding to different radii in the dual geometry Thus, with 
the holographic Wilsonian renormalization group (RG) approach [15, IE, [T7|], it is also 
interesting to study the holographic hydrodynamics at a finite cutoff surface directly 



14, 18 



It was proposed in ]21| that there is a mathematically precise relationship between 
the unforced incompressible NS equations in p + 1 dimensions and vacuum Einstein 
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equations in p+2 dimensions and their solutions. The dual geometry has an intrinsically 
flat timelike cutoff surface whose extrinsic curvature is identified as the stress energy 
tensor of the dual fluid. This relationship has been further developed in the literature 
such as [22, 23, |24|, |25|, p7[ . In [^TJ, a gravitational shock wave was introduced to 
stir the fluid and then left to evolve according to the unforced NS equations. In order 
to study the stationary NS turbulence, it is better to introduce the external random 
source fields 0, |28| . For example, a dilaton field was added to the bulk gravity in fll2"|. 
In this case, a perturbed gravity solution with a slowly varying dilaton leads to a slowly 
varying force term in NS equations. Another example is [[13[], where the force terms 
come from the slow variation of the boundary background in the holographic context. 
It turns out that a simple forced steady state shear solution to the forced NS equations 
becomes unstable and may translate into turbulence at high enough Reynolds number. 

In the fluid/gravity correspondence, the derivative expansion method [1C] leads 
to the equations of motion of relativistic fluid, and the equations reduce to the NS 
equations in the non-relativistic limit [ |T3| , p9 . Actually, the NS equations can also 
be obtained via taking the non-relativistic expansion method directly ]2T], T2\. This 
method can be used not only for dual fluid at the AdS boundary, but also at a finite 
cutoff surface with flat induced metric []18| . In the latter case, the bulk geometry is not 
required to be asymptotically AdS. While in the asymptotically AdS case, the perturbed 
gravity solutions with Dirichlet condition at the cutoff surface can be mapped to the 
perturbed gravity solutions without the cutoff surface pO, pll. To study the holographic 



fluid with external forces in non-relativistic limit, instead of using the induced metric 
perturbations method proposed in O, RIL we keep the induced metric flat and add 



external matter fields in the bulk as the source terms to the dual fluids. 

It was shown in |TJ| that the radial component of Einstein equations is equivalent 
to the isentropy equation of the dual fluid. Using a more generic static metric, we 
show that this statement is also true in the Gauss-Bonnet gravity. Using the non- 
relativistic expansion method, we give the procedure to obtain the perturbed solutions 
of the bulk gravity with a finite cutoff surface up to the second order of non-relativistic 
expansion parameter, and corresponding NS equations of the dual fluid at the cutoff 
surface. It turns out that the shear viscosity over the entropy density of the fluid dual 
to the Gauss-Bonnet gravity does not run with the cutoff surface. This part acts as the 
service to introduce the non-relativistic expansion method. In this paper we mainly 
focus on an Einstein-Maxwell-dilaton system with a negative cosmological constant and 
pay attention to the external force terms in the dual non-relativistic fluid. The external 
force terms come from the Maxwell field and dilaton field in the system. Note that the 
forced fluids at the AdS boundary have been discussed in the Einstein-dilaton system 
12| and Einstein-Maxwell system |32], [33], |3~4| , respectively. The perturbed solutions 
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have been obtained to the second order of the derivative expansion. We consider the 
Einstein-Maxwell-dilaton system and obtain the perturbed solutions up to the second 
order of the non-relativistic expansion with/without the cutoff surface. Associated 
forced NS equations are also derived in both cases. The concrete expressions of external 
forces of the dual fluid are given. The results show that the Reynolds number of the 
dual fluid becomes larger and larger when the cutoff surface approaches the horizon of 
the background black branes. 

This paper is organized as follows. In Sec.||, we start with a generic static black 
brane metric and obtain the perturbed solutions with a finite cutoff surface in the 
Gauss-Bonnet gravity by using the non-relativistic expansion method. This section 
acts as to fix the notations in this paper and to introduce the non-relativistic expansion 
method. In Sec.[3], we consider the Einstein-Maxwell-dilaton system with a negative 
cosmological constant and discuss the dynamics of dual fluid on the AdS boundary 
in non-relativistic limit. We generalize the discussions to the case with a finite cutoff 
surface in Sec.[|. The conclusions are given in Sec.||. 



2. Holographic fluid at a finite cutoff surface 



This section is a generalization of discussions in |L8[ on the thermodynamics and hy- 
drodynamics of dual fluid at a finite cutoff surface. Slight differently, we start with a 
more general static metric and work in the intrinsic coordinates on the cutoff surface 
directly. 

2.1 Thermodynamics of the dual fluid at the cutoff surface 

To study the fluid in a (p + l)-dimensional flat spacetime, we consider the generic 
(p + 2)-dimensional static black brane background 

ds 2 p+2 = -g tt (r)dt 2 + g„{r)dr 2 + g xx {r)bijdx % dx 3 , = 1,2, ...,p, (2.1) 

where the metric components are functions of radial coordinate r only. We assume the 
metric has a well defined event horizon at r = r^, where gu{f) has a first order zero 
9tt( r h) = 0, and g rr (r) has a first order pole g^^h) — HI5fl - For example, the ingoing- 



Rindler form of flat spacetime and the black p-brane solutions in asymptotically AdS 
spacetime have the form [jTJJ . Using the Eddington-Finklestein coordinate r defined by 
dr = dt + a/ g rr ( r ) 1 9tt( r )dr |y| , we can rewrite the metric (]2.1|) as 



ds 2 2 = 2yJ g tt (r)g rr (r)dTdr - g tt (r)dT 2 + g^^S^dx'dx 3 , (2.2) 
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which has the translational invariance in r and x % directions. We can always introduce 
a finite cutoff surface E c at r = r c outside the horizon with the intrinsic coordinates, 

x a ~ (f , x l ), as 



x° 



{a,b...} = 0,l,2,...,p. 



= r = y/gtt(r c ) r, x l = V 'g xx (r c ) x\ 
Then the associated bulk metric ( |2.2| ) becomes 

ds 2 p+2 = g r r(r)dr 2 + lab (r) [dx a + N(r)5 a f dr] [dx b + Af(r)5 b dr] 
where ^ a b{ r ) and N(r) are given by 



j ab (r)dx a dx b 



y/9tt(r c )9rr(r) 



gtt{r c ) 9xx{r c ) y/gu(r) 

And the induced metric with intrinsic coordinates on E c is simply given as 

1ab( r c)dx a dx b = fj ab dx a dx b = —df 2 + S^dxtdx, 3 . 



d~s 2 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



The Brown- York stress energy tensor Y evaluated on the cutoff hyper surface S c 
is proposed as the stress energy tensor of the dual fluid pif . It has a close relation 
with the extrinsic curvature tensor K ab = \C^ ab {r)\ r=rc of E c , where Cfj is the Lie 
derivative along the unit normal N A of the hypersurface. To study the thermodynamic 
properties of the dual fluid at the cutoff surface, we begin with the re-scaled metric 
( |2.4| ) with a Killing horizon located at r = r^. The local temperature T (r c ) on S c , 
which is identified as the temperature of dual fluid []14 |, meets the Tolman relation with 
Hawking temperature Tjj of the black brane metric 



T (r c ) 



T, 



H 



Tn = h m 



gtA r ) 



(2.7) 



V9tt{r c ) r -*r h 4my/g tt {r)g„(r) 

To discuss the local entropy density of the dual fluid, we consider a quotient of the 

x % + £on l , with a characteristic length 



x 



general geometry (|2.2|) under shift of x % |14 
£o and ri 1 G Z. Equivalently, using metric (|2~4]) , the spatial BP on E c turns out to be a 
p-tours T p with r c -dependent volume 



V P (r c ) 



g p J x 2 (r c )v 



Vo = f , x'^x l + £ y/g xx {r e )ri 



(2- 



As the total Bekenstein-Hawking entropy 5* is fixed, we can identify the dual fluid's 
entropy density as 



S 



1 



gxx{r h ) 



V P (r c ) 



S 



V P (r h ) 



AG 



p+2 



(2.9) 
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where S is the Bekenstein-Hawking entropy of the black brane. 

As a calculation example, we consider the Einstein-Hilbert action with a non- 
positive cosmological constant A < 0, as well as the Gauss-Bonnet term £gb = R 2 — 
ARabR AB + RabcdR ABCD , and the contribution from matters Cm, 



1fi \ [ d p+2 x v ^(R-2A + aC GB + C M ) 
167rGp +2 J 



(2.10) 



where a is the Gauss-Bonnet coefficient with the same dimension as square of length. 
The Gauss-Bonnet term is a topological invariant when p = 2, and will modify Einstein 
gravitational field equations when p > 3 [f}5|, [36| . We will generally demand p > 2 and 
set 16nG p+ 2 = 1 throughout this paper 1 . Varying the action fl2.10|) with respect to 
the induced metric •y ab with appropriate surface terms |37|, [39], f|0|, we can get the 
Brown- York stress energy tensor Y evaluated on the cutoff hypersurface S c of the 
Gauss-Bonnet gravity. Because of the flatness of the cutoff surface, its final form is 
given by 



77, 



BY 
ab 



Kr] ab - K ab - 2a 3J ab - Jr] ab ) + Ci], 



I a I, 



with 



J ab = - { 2KK ac K* b + KdK Kab - ™ ac K*K & - K 2 K ab) 



(2.11) 



(2.12) 



where K is the trace of extrinsic curvature tensor K ab of S c , and C is an ambiguous 
constant. Substituting the metric (|2.4j ) in the Brown- York stress tensor (|2.11|) , we can 
get the stress tensor of the holographic fluid dual to the Gauss-Bonnet gravity as 



T a BY dx a dx b 



Pg\Tc) 



2C 



df 2 + 



p G (r c ) + 2C 



(2.13) 



with constant pressure p G (r c ) + 2C and energy density p G (r c ) — 2C, where 



P G (r c ) = G p (r c )p E (r c ), 



P 9 xx {r c j 



V9rr{r c ) 9xx(r c 



1 9L(r c )\ (2-14) 



y/g rr (r c ) \9tt(r c ) 9xx(r c ) 



The sum u G (r c ) is independent of the constant C. The Gauss-Bonnet term corrections, 
compared to Einstein gravity with \g = {p ~ l)(p — 2) a, are given by 



Q p {r c ) 



A 



G 



9'xx(r c ) 
69rr{r c ) \9xx(r c ) 



A 



G 



9'xx{rc) 
^9rr{r c ) \9xx(r c ) 



(2.15) 



1 Most of the results in this paper have been checked by Mathematica up to p = 5. 
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Using the above results, in general we have the thermodynamic relation 

T (r c )s (r c ) = u G {r c ) - ty,(r c ), ty,(r c ) = ^ ^ (r c )n^(r c ), ^.16) 



where m = 1,2,..., and every n^j(r c ) corresponds to different kinds of charge density 
with corresponding conjugate chemical potential /im (r c ). The total entropy is given by 



S = s (r c )V p (r c ) - V " lr ' ] 



T (r c ) 



d r S = 0. 



(2.17) 



This isentropic equation d rc S = can be considered as either an adiabatic thermo- 
dynamic process of the dual fluid or a holographic renormalization group flow [H|. 
Varying the action ( p.!0|) with respect to metric g AB , we can obtain the equations of 
motion for the Gauss-Bonnet gravity in the re-scaled p+2-dimensional bulk coordinates 



W AB {r) = E AB {r) + aH AB (r) - - f AB {r) = 0, 
with the definitions 



{A,B...} = {r,f,&} (2.18) 



1 5j0 

E AB (r) ee R AB - -R g AB + Ag AB , f AB (r) = -2j^ + g AB C M , 



(2.19) 



H AB (r) = 2(R AC deR b —2R ACB dR —2R ac R b + RR AB ) — -g AB £>GB, 

where T AB [r) is the stress energy tensor of bulk matters. If we assume the metric 
fl2.4j) solves equations ( [2.18|) , the non-zero components of the stress energy tensor are 
f rr (r),fff(r),f ss (r) = f&&(r), and f fr (r) = f rf (r) = N(r)f ff (r), respectively. Fol- 
lowing the procedure in jT3| and using metric (2.4j), we can introduce an arbitrary null 
vector ( tangent to the hypersurface r = r c with time component d? , such as 



( A d A = d f - % g AB ( A ( B = 0, ( A N A = 0, 

with N the unit normal of the cutoff surface. It could be checked that 

V p (r c ) 



Or 



T o( r c) 



T, 



H 



■V9rr(r c )g tt (r c ) 2 ( A ( B 



E AB (r c ) + aH AB (r c) 



(2.20) 



(2.21) 



In general, the Cl (r c ) term would satisfy the following relation 2 



* \'grAr,)g !l {r. ) ( A ( B f AB (r c 



T, 



H 



(2.22) 



2 For an example see Sec.0 of this paper or |L4| 
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Using (|2.17|) , (|2.18|) and the above two equations, we have 



d rc S = Y^^g rr (r c )g tt {r c ) 2 ( A ( B W AB (r c ). (2.23) 

J-H 



As ^/grr(rc)gtt(r c ) } V p (r c ) and T# have neither zero nor pole even when r c = r^, we 
arrive at 

d rc S = « C A C B ^AB(r c ) = 0, (2.24) 

which implies the equivalence between the isentropy of the RG flow and a radial grav- 
itational field equation of the Gauss-Bonnet gravity. 

2.2 Incompressible Navier-Stokes equations from gravity 

Keeping the intrinsic metric ( |2.6| ) of S c flat, we can take two linear diffeomorphisms 
based on metric ( |2.4|) : a transformation of the radial r and a Lorentz boost in (f, x l ) 
coordinates. The first is the transformation of r and the associated re-scaling of (f, x l ) 

as 

. ; / \ , ~ / 9ttyc) ~i ~i \ Qxxyc) / n nz\ 

r^k(r), r-»r, , x l -> x\ ftt yp 2.25 

where fc(r) is a linear function of r, whose form would be chosen via the global symmetry 
of geometry Q2.4|) . The metric ( p.4|) is then transformed into 



ds* _> ds 2 = 2 ^ 9t f )9r ^ ) dfdf - ^Prdf 2 + to^Sydx'dxi. (2.26) 

V9tt(r c ) gtt{r c ) 9xx{r c ) 

where we have introduced the re-scaled coordinate r = k(r) and the notation f c = k(r c ). 
We will work in the (r, x a ) coordinates directly in this paper. 

The other diffeomorphism is the Lorentz boost with a constant boost parameter fa 
in the x a = (f , x % ) coordinates 

f _> -u a x a , x l -> h\x\ n« = (- 7j 8«, 5) + ( 7 - 1)^) , (2.27) 

where we also have defined the (p + l)-velocity 

u a = 7s (-l, A), 7, = (1 - 2 )~* , /3 2 = PiF = ^-W- (2-28) 
As these parameters are all constants, the associated transformations can be expressed 



as 



df -»■ -u a dx a , dx % -> n J a cfe a , S^d^dSP -> 5 ij n l a n{dx a dx b = V ab dx a dx b , (2.29) 



-7- 



with the projection operator V a b = Vab + u a Ub = bijn % a n\. Then the metric (|2.26|) 
becomes 

df p+2 ds 2 p+2 = g rr {r)dr 2 + % b {r) dx a + N{f)u a df dx b + M(?)u b df , (2.30) 
where we have used the ADM-like decomposition at the constant f surface, 



~ 9tt(r) ~ ~ , 9xx{f) y Kn~\ V9tt(r c )grr(r) 

lab{r) = -prrUaUb H j^-Vab- N (r) = — - . = = — . (2.31) 

9tt(r c ) 9xxir c ) \/gjJ) 

After the two diffeomorphisms, the metric ( j2.4|) and ( 2.30Q still solve the same gravity 
field equations. The cutoff surface r = r c is equivalent to f = f c , we will firstly work 
with the new radial variable r , and then transform back to r via r = fc _1 (f). 

In general, the Brown- York stress energy tensor on the cutoff surface r = r c 
corresponding to metric (|2.30|) will turn out to be taken the following form 

7^t y (^ c ) = p(r c )u a u b + P(r c )V ab , u(r c ) = p(r c ) + p(r c ), (2.32) 

which could be identified as the stress tensor of an ideal relativistic fluid j^ dea Q with 
(p + l)-velocity u a in flat space-time. In the Gauss-Bonnet gravity case ( J2.ll ), the 
cutoff dependent energy density p(r c ) = p G (r c ) — 2C and pressure p(r c ) = p G (r c ) + 2C 1 . 
In general here C is an unfixed constant, but to obtain a finite result when the cutoff 
surface goes to the AdS boundary, C can be fixed fll8| . 

Let us pause to recall the hydrodynamical description of microscopic field dynam- 
ics in flat spacetime. It applies when the correlation length of the fluid l cor is much 



smaller than the characteristic scale L of variations of the macroscopic fields [41 



Via dimensional analysis, l cor ~ 1/T C and 1/L ~ d^, where T c is the characteristic 
temperature of the fluid and would act as the coordinates dependent parame- 
ters, such as T{x) and u a (x) 3 . One introduces the dimensionless Knudsen number 
Kn = lcor/L ~ ^d a ~ e 1 to expand the stress energy tensor of relativistic fluid in 
flat background fj ab , 

oo 

T ab (x) = ^tf^), t} b 0) {x) = u{x)u a {x)u b {x) +p{x)f) ab , (2.33) 

n=0 

and T^ b n \x) ~ (Kn) n . To take the non-relativistic limit as in jlTJ, we recover the speed 
of light c in (|2.33|) and introduce the normalized pressure P(x) as 

u a (x) = lp {x) (-1, 0i(x)/c) , l0 (x) = (1 - (3 2 (x)/c 2 y\ P 2 (x) = k(x)P(x\ 
u a {x)d a = 'j{x)— + ^{x) d,, 



c ' tpK ' c % ' T(x) u(x) p(x) p(x) c 2 

(2.34) 



3 We have used (x) to denote the function arguments (x a ), and would use d a to denote dx 
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When c — > oo, the energy momentum conservation equations of the ideal fluid, d a T^ (x) 
0, leads to the non-relativistic incompressible Euler's equations, 

diP(x) + drfiiix) + ^{x)djfii{x) = 0, diP\x) = 0. (2.35) 

Instead of the c — > oo limit, if we assume the small velocity parameter to take the same 
limit as the Knudsen number, i.e. /3/c ~ Kn ~ e, it would be equivalent to set c = 1 
in ( |2.34| ) with the following scalings 

di~e, B T ~e\ Pi(x)~e, P{x)~e 2 . (2.36) 



These non-relativistic scalings are named as the BMW limit [13, M. To obtain the 



forced NS equations, the dissipative part of the stress energy tensor (|2.33|) is required 

T ab (x) = t^(x) + t} b diSS \x), t^ ss \x) = f t f(x) + 7%\x) + ... . (2.37) 

For example, in the Landau frame u a (x)T^ b ss (x) = 0, the first order dissipative com- 
ponents could be written as 

7^6 , (£) = -2fj(x)a ab (x) - ((x)6(x)V ab (x), V ab {x) = fj ab + u a (x)u b (x), 

P h (x)~ ~ ,~ ( 2 - 38 ) 

Zab(Z) = V™{x)V%{x)d {m u n) {x) - ^-Le(x), 9{x) = fj ab d a u b (x), 

p 

where fj(x) is the kinetic shear viscosity and ((x) is the bulk viscosity which will vanish 
for conformal fluids. Usually, they behave the same as local temperature T(x) in ( j2.34|) . 



The first order dissipative hydrodynamics satisfies the dynamical equations d a [T} b \x) + 
7^(5;)] = fb{x) if some external source term appears. In the non-relativistic limit 
(|2.36|) , if we further assume fi(x) ~ e 3 , ff (x) ~ e 4 , it would lead to the non-relativistic 



forced incompressible NS equations at order e 3 , 

BiP(x) + drPiix) + Fiftd^x) + vffl&djPiix) = ff ] (x), d^\x) = 0, (2.39) 

where v(x) = r](x)/£d(x) is the dynamical shear viscosity and f} u \x) = fi(x)/Q(x). 
The NS equations have the scaling symmetry ( |2.36| ), as shown in [|13|, 



Next we will derive the incompressible NS equations from the gravity side. To 
get the dissipative part of the dual fluids, we need to perturb the geometry ( |2.30| ) by 



using either the linear response method (eg. ||15| ) or the perturbative expansion method 
developed in flOf , where a perturbative procedure to solve Einstein's equations order 
by order in the boundary derivative expansion was proposed. In the case with a finite 



cutoff surface, we also can employ the procedure under the non-relativistic limit [22, 18 . 
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Regarding the transformation parameters in linear function r = k(r) and the p- velocity 
(3 l in the boost transformation as functions of the hypersurface coordinates (f,x l ), we 
can solve the gravitational equations order by order in the non-relativistic perturbative 
expansion. Based on ( |2.36| ), these transformation parameters are also regarded as small 
quantities with appropriate scaling symmetry as 



d r ~ e°, d % ~ Pi(x) ~ e 1 , d T ~ 5fc(r, z) ~ ~ e 2 



(2.40) 



where £fc(r) = f — r ~ e 2 would provide the pressure perturbation. Using the formal 
Taylor expansion such as g tt (f) = g u (r)+g' tt (r)6k(r), and define g fr (r) = \/gtt(r)g r r(r)/g t t(r c ), 
we can consider the metric fl2.30| ) with coordinate dependent parameters as ds 2 ^, and 
expand it up to order e 2 , 4 



d s (o) — + 2gfr( r )drdr + 



gtt{r) 



- 2g fr (r)/3i(x)dx l dr + 2 



gu{r c ) g X x{r, 

9tt{r) g xx (r) 



(5^ j doc dec 



gtt{r c ) g xx [r, 
9tt(r) g xx (r] 



(3i(x)dx l dr 



_gtt(r c ) gxx(r c ) 



+ gfr(r)^flsk(r)dfdr + 

grrir) 

+ 2g fr {r) [jfe'(r) - l]dfdr + 
+ 0(e 3 ), 



9 ^ r) 5k(r) - 9 -4i8k{r c 



9tt[r) 
9xx(r) 

9xx(j"c, 



9tt{r c 



g'xx( r ) 5k ( r ) 
gxxir) 



9xx\Xc 



gf r (r)dfdr — ® tt \ \ dr 2 
9tt{r c ) 

-5k(r c ) dxidx\ 

(2.41) 



where Sk(r) and /3j(x) and are all functions of x a now, we will omit this notation (x) 
henceforth. The coordinates dependent metric (|2.41|) is no longer a diffeomorphism of 
metric (|2.30|) . Under the perturbative expansion, it turns out that (|2.41 ) only solves 
the same equations of motion of gravity up to e 1 . To solve the equations of motion 
up to e 2 , some constraint equations and new correction terms to the bulk metric are 
needed at this order. With the same way, the equations of motion of the system can 
be solved order by order in the non-relativistic expansion parameter e ||22|| . And the 



corresponding Brown- York stress energy tensor which is identified as the dual fluid's 
stress energy tenor, can also be obtained at the desired order. In this paper we solve 
the equations of motion up to the order e 2 . It turns out that the non-dissipative part 



in this paper, we use the scripts (0), (1), ... to denote the order in the derivative expansion, and 
scripts (e), (e 2 )--. to denote the order of the non-relativistic expansion parameter e. 
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is still given by ( |2.32|) in the non-relativistic limit. Up to e 2 , we have 

[p(r c ) + Cj Q (r c )(3 2 ] df 2 - 2Q (r c )(3 i dx i df 



t^dx a dx b 



+ [p(0<% + %( r c)A/3j] dx l dx j + 0(e 



(2.42) 



where u (r c ) = p (r c ) + p (r c ), and 



P( r c) = P (r c ) = p (r c )+p' (r c )8k{r c ), p(r c ) = p (r c ) = p (r c ) + p' (r c )5k(r c ). (2.43) 

The equations of motion of the bulk matters should be solved in the same procedure. 
Once we get the solutions of both gravity and matters up to order e 2 , the constraint 
equations of gravity at order e 3 are just the forced incompressible NS equations (|2.39|) 
of the dual fluid. 

As a calculation example to perturb the geometry ( |2.30| ), we again assume the met- 
ric solves the equations ( |2.18| ) of Gauss-Bonnet gravity with a non-positive cosmology 
constant A. The metric Q2.41Q leads to W^ N = 0, and the constraint equation at order 
e 2 turns out to be 



-8^ = ^)8^ = 0, 



(2.44) 



where £> (r. 



= uj g (r c ) is nonzero outside the horizon. Thus it leads to the incompress- 
ible condition d^ 1 = of the dual fluid at the cutoff surface. 

To solve gravitational field equations at order e 2 , we need to add corrections to the 
metric (|2.41|) . It was shown in [10] that due to the spatial SO(p) rotation symmetry of 
the black brane background, one has the decoupled equations of SO{p) scalar, vector 
and traceless tensor perturbations. For example, if we turn off the tensor perturbations 
of the matter stress energy tensor Tab{^), the tensor corrections to the metric, at order 



arc 



ds 2 {1) 



2JP(r) 



9a 



dij dx l dx J + 



<?ij = d(ifo) - -5ijd k (3 ' 



Then the traceless tensor parts of the gravitational field equations 
second-order ordinary differential equation 



(2.45) 
lead to a 



i + 

V9rr(r) 



Ga(r) 





l 




lr \/gtt{y) 







dy, (2.46) 



where two integration constants have been fixed and 



/ 



l-ot(p-2) 



9tt(r)g. 



(p-3)/2, 



9 xx {r) 



V 



[gtt{r)grr{r)}gxx 1)/2 (r 



(2.47) 
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In (|2.46|) , one integration constant is the upper limit of the integration, which is chosen 
to be r c via the Dirichlet boundary condition at the cutoff surface. The other integration 
constant is chosen to be Q g (rh) to cancel out the first order zero in the denominator 
V ' 9tt(jh)/ a/ g r r(fh) = of the integrand in ( |2.46| ), where 



Gg(r h ) 



lim Q g 

r^r h 



9 P J x 2 (r h ) 



^p-2)ixT H g' xx (r h ) 
\/gtt(r h )g rr (r h ) 9xx (r h ) ' 



(2.48) 



and a/ gtt(rh)grr( r h) is a finite constant due to our definition. The Hawking temperature 
Th is given in Q2.7|) , so that A a is a constant determined by the metric ( p.l|) at the 
horizon. 

With the perturbed metric ds 2 ^ + ds 2 ^, corresponding Brown- York tensor becomes 



tKo) , ^(i) 

'ab 1 " 



ab ' 



where the symmetric traceless components of T2p are 



r {1) 

>ij 



-2 



9xx(r c 



1 + y^P(r c 

V9rr\r c ) 



_ G g (r h ) 
9xx (r c ) 



-2fj(r c )a i:i . (2.49) 



Here fj (r c ) is defined as the kinetic shear viscosity of fluid dual to the Gauss-Bonnet 



gravity with the geometry ( [2.1|) . Using the entropy density s (r c ) given in ( |2.9| ) of the 
fluid at the cutoff surface, we have 



a p/2 (r 
gxx yi c 



fj(r c 



-u 

4tt 1 



aA r 



(2.50) 



One can see that the shear viscosity over entropy density of dual fluid does not run with 
the cutoff surface. This match with the results in jn| and [1~2"|, where a hypersurface 
in the given solution is introduced. Here we impose the Dirichlet boundary conditions 
IfLil, |i~8|1 and find the corresponding perturbed bulk solutions with a flat induced hy- 
persurface. Using the formula ( [2.50|) , one can easily obtain the ratio of shear viscosity 
over entropy density for some fluid once the dual gravity background is known. Take 
an example, using the charged AdS Gauss-Bonnet black brane background solutions, 
one has |19| 



9xx(r) 



2 9 

p, 9tt(r) = g-\r) = — 



4A 



G 



fj(r c ) 
s(r c ) 



1 

47T 



1 - 2 



G 

\<; f p+1 
p — 1 



A G = (p-l)(p-2) Qtj 



(2.51) 



where a is the Gauss-Bonnet coefficient, and £ is the radius of AdS spacetime. When 
t — > oo, the negative cosmological constant A = — p(p + l)/2£ 2 — > 0, the black brane 
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solution will degenerate into the ingoing Rindler space via some coordinate redefini- 
tion fi\M, and the Gauss-Bonnet term will not contribute to the shear viscosity any 



more. This has been shown via using the ingoing Rindler metric of flat spacetime [ I3| . 
Here our formula in (|2.48|) lead to this result directly because g X xij') = 1 =>■ A a = . 



3. Holographic charged fluid at AdS boundary 

It is interesting to give the concrete expressions of some external forces of holographic 
fluid. For this, let us consider the Einstein-Maxwell action with a dilaton field $ and 
appropriate boundary terms and counter terms as follows. 

where the negative cosmological constant A = —p(p + l)/2£ 2 , the Maxwell field Fmn = 
2V [M -4jv], and F 2 = F MN F MN , (V$) 2 = V A /$V M $ , where {M, N...} run over the 
bulk coordinates. In what follows, 167rG p+2 = 1 will be adopted. From the action we 
can get the equations of motion of the Maxwell field and the dilaton field 

W N = V M F% = 0, =V 2 $ = 0, (3.2) 

as well as the gravitational field equations 



D 9m N 
tf-MN ;r~ 



R + Afl'Miv - - (Tmn + t mn) 



T ( A ) - TP TP p 9MN ^2 T ($) _V7 ^ a 3MN_ (X7 ^2 



(3.3) 



from which one has 



W MN = Run + ~~~~jpr~9MN - ld M $d N $ - ~ ( F MP F/ - ^-g MN F' 2 ) = 0. (3.4) 



2p 

The system has a class of exact solutions with p + 2 parameters as 



ds = -lu^dx^dr + — [V ' (r)u^u v dx^dx u + V^] dx^dx 1 
A = -U(r)uudx>*, $ = $ - 



(3.5) 



where {//, v, ...} run over the boundary coordinates. This is a charged boosted black 
brane solution in Eddington-Finkelstein coordinates with a constant dilaton <3> , and a 
constant unit normalized velocity w M , 



u, 



7„(-l, Vi), 7, = (1 - v 2 ) 5 , v 2 = vy = SijvW. (3.6) 



13 



and rf^u^Uu — — 1. In the metric of the solutions 

m£ 2 



V{r) 



+ 



2 1)2 



yP+1 iy2p 

Here the notations in f03l have been used, 



Q 



P+l 2P 
1 ^ V r 2p • 



(3.7) 



V(r h ) = 1-M + Q 2 = 0, 



M 



m£ 2 
3+i ' 



Q 



q£ 



(3.8) 



where r = is the horizon location of the charged black brane solution. The Hawking 
temperature of the black brane is given by 



r 2 V(r) 
~1 2 ~ 



[(p + l)M - 2p Q 2 ] 



r=r h 



4vr£ 2 ' 



(3.9) 



In the gauge field part, we have chosen the gauge that A r = and 



U{r) 



Mo 



2p 



P 



(3.10) 



where 6oMo is the value of the gauge potential at the AdS boundary. Usually bo = 1 
is required to get a well defined gauge field at the horizon §441 [45|, but this is not 



necessary in the perturbative fluid/gravity correspondence |32 
3.1 Charged fluid in non-relativistic limit 



To perturb the metric (|3.5|) , we take the same procedure via regarding the parameters 
in the metric and gauge field as functions of boundary coordinates —> u^(x), r h — >■ 
r^(l + P{x)). 5 Without loss of generality, we take the associated replacements as 
m — > m (1 + P(x)) p+l and q — > q{\ + P(x)) p to keep M and Q as two constants. 
The scalar field is replaced as $ — > 0(x M ). Then the solutions ( |3.5|) will not solve the 
equations of motion Q3.2|) and ( |3.4j) any more. Using the so called BMW limit at the 
boundary [I 



21 



8 r 



di ~ t> j ~ (9j0 ~ e , <9 T ~ P ~ e 2 



(3.11) 



where we have added an assumption that di<p ~ e , we can solve the equations of 



motion order by order in the small parameter e [£2|, [L8| . Actually the three parameters 
(m, q, r^) relate to each other via U(r^) = 0, so that two of them are independent 



5 we briefly denote the function variables (x^) as (x) and some of them in the following equations 
will be ignored. 
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parameters, for example, one can take q and as two independent parameters. We 
focus on the forced NS equations in this paper, so 5m ~ 5q ~ 5r h ~ e 2 have been 
assumed in the non-relativistic limit. In the boundary derivative expansion of the 
relativistic fluid, the solutions of the Einstein-Dilaton system [fL2], and the Einstein- 



Maxwell system []32], have been obtained up to the second of the derivative expansion 
parameter. We can extract the relevant terms up to the second order in our non- 
relativistic expansion parameter e. To compare with their results, we write our the 
non-relativistic perturbative solutions up to e 2 in a covariation form, 

ds 2 = - 2u IM (x)dx> l dr + — [-V (f) u^{x)u u {x) + V^{x)\ dx p dx u 



+ 



„ + 2lT{r)a ilv + £ 2 J r (r)(^)^] dx»dx h ', (3.12) 



A =A^(r,x)dx» = [Af{r,x] + A^\x)] dx», 
$ =0( x ) + tT(j) (u»d^) + t 2 ^{r){d»d^). 

Although higher order terms such as u T Ui ~ e 3 are reserved, they would not contribute 
to the calculation up to order e 2 . In (|3.12| ), V^ u {x) = rj^+u^xju^x), r = [1 — P(x)] r, 



1(9)/ 



and Alt" [r,x) 

an extra electromagnetic field, which would provide extra forces and F f 



[1 + P(x)] U(r)Ufj, is the perturbed and boosted potential. is 



(ex) 



are chosen to meet with the appropriate non-relativistic expansion in |p~4 

^3 



equations of motion (|3 
provided 

T{r) 



e , as well as the velocity u. 



and 



(— 1 — v 2 /2, Vi). This solution solves the 
4|) up to order e 2 , if the following expressions are 

(y p - <) 



y p+2 v( y ) 



dy, 



a 



fJtV 



-p v 

d^u v) —6, 



P 



where 9 = rj^d^Uv = V^d^u^, for the dilaton shear perturbations 



p — 1 



(v 



p-i _ _p-i 



y p+2 V(y) 



dy, 



\<y<t 



I /IV 



V, 



V 



and for the dilaton scalar perturbations 

/C(r) 



(3.13) 



(3.14) 



(3.15) 



2p(p — l)r 2 ' 

We can generalize the dual boundary fluid in [O, B3J to the (p + 1) dimensional case 
with p > 2, where the stress tensor 7j" and current are given by 



r p+l 
lim - 



2p. 



2C 



2(K a ph e *h* - K$) - -fK + T G " 



r 



P+i 



lim 



NrF^, 



p — 1 

(G 



p — 1 



(3.16) 
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where N is the outward pointing unit normal of the regulated boundary, D is the 
covariant derivative associated with the boundary metric and the corresponding 
Einstein tensor G%, as well as the effective counter tensor (G$)% of the dilaton field 
pT7| , Substituting our solutions ( |3.12[ ) in ( |3.16| ), the Brown- York stress tensor 



and the induced current at the AdS boundary are given by 



M h 



r 



p 



2 [1 + (p + 1)P] [(p + 1) u fJ 'u u + rT\ ~ 2 ^ a"" 



0"„ 



J' 



pur 



n 



(?) 



1 + pP)n , 



n 



-C/'(r) 



P l " (3.17) 

_ q^2p(p- 1) 



where n.( 3 ) is the induced charge density. Note that here (cx^)^ is included in the 
stress tensor (p. 17 ). These terms relating to the dilaton field were considered as the 
dissipative parts of the dual fluids's stress tensor in [0, or the perturbations in the 



membrane paradigm JT5|. In this paper, however, we focus on the forced NS equations 
and move them to the right hand side of the constraint equations as additional sources, 
similar to what has been done for the boundary metric perturbations ]13|, j3l| . Then 
we can redefine the stress tensor of the dual fluids with first order dissipative term 

^,p+ 1 



7>> = co [1 + (p + 1)P] 



u^u u H — 

p+1 



2i] a 



UJQ 



(p + l)M-£^, (3.18) 



where r] Q is the holographic shear viscosity of the dual fluid. With the corresponding 
entropy density Sq, we have 



so 



Vo 
so 



1 

4tt' 



(3.19) 



4G p+2 Ip 

and the thermodynamic relation ThSq = ooq — no/io still holds. Furthermore, we can re- 
define the following stress tensor by moving out the terms associated with the boundary 
chemical potential 



7$ = T H s [1 + (p + 1)P] 
This form will be used later. 



U^Uv + 



p + 1 



- 2 770 <r, 



(3.20) 



3.2 Forced Navier-Stokes equations 

Using Gauss-Codazzi-Mainardi relations near the boundary [[L2], we can obtain the 
constraint equations of gravity 

2D A (K CD h CD hi - Kg) = -2R CD h c B N D = -D B §V D <b N D - F CP F p D h% N D , 

(3-21) 
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where Kab, N a and Kab are respectively the induced metric, the outward pointing unit 
normal and the extrinsic curvature of the constant r hypersurface, with the covariant 
derivative D. When r — > oo, transforming these results into the boundary metric r\^ v 
of the dual fluid, we have 



r P+i 



8 T M 

U H 1 V 



lim 



N C V C $ + 



Introduce the gauge field tensor Fjtv = 2d[ fl A^ which comes from the boundary 



p — 1 



D*Q ) A,$ + N c F CD F y 



D 



(3.22) 



chemical potential A 



(bd) 

dbg] 



-b p (l + P)u^ , and denote the background gauge field 

lOs) _ a( M ) i A( ex ) 



tenor as F^f = 2d\ tl A u ^ , where A^ = A^ + A^ , the (p+ 1) constraint equations 
up to e 3 become 



-d v 



-£ (u»d a cb) + 



p — 1 



p-1 



^i^, (3.23) 



where F^t = F^ x> + F^™' , and the constraint equation of Maxwell field reads d^J^ = 
n^d^u^ = 0. Using the redefined stress energy tensor ( |3.18|) , we get the forced NS 
equations at order e 3 , 



dnS ) = + fi q \ fi q) = -n iq) u^ 



where is just the Lorentz force of the charged fluid and fy 9 ' is the external force 
from the dilaton field up to e 3 



(4) 



(3.24) 



o p— 1 

P r/l -. (3.25) 



p{p - I) 



The Lorentz force up to e 3 can also be divided in to two parts as 



f ( q) ^ f (e X ) + / CW) j = _ mu *Fjfi*) 7 f V-, = _ n{q)U r . 

In addition, we notice that the following relation holds at order e 3 

P[T$ ) -T$]b = fS*> ^ 8n$ = f<f> + (3.27) 

Thus, if we define the fluid stress energy tensor as (|3.20|) , with 6 = 1, we can see that 
the external forces only come from the dilaton field and the external electromagnetic 
field. It would be more clear to see the non-relativistic results by neglecting higher 



l Fl d) . (3.26) 
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(3.28) 



order terms of the solutions in the covariant form. The solutions in ( [3.12|) become 

ds 2 = + 2drdr + — [-V (r) dr 2 + tfyCteW] - 2v i dx i dr - 2— [1 - V(r)] v i dx i dr 
+ v 2 drdr + T — [rV'(r)Pdr 2 + (1 - V (r)) (v 2 dr 2 + v i v j dx i dx i )] 

2 

+ T - [fJC^dr 2 + {2lF{r)<Tij + ^(r)(t7 ) y ) dx'dx?] , 
A = [U(r) (1 + P + v 2 /2) - rPU'(r) + A^ x) ] dr + \-U(r) Vi + A { - x) 
{ $ =0( x ) + tF{r) (<9 r + v%(£) + £ 2 ^(r)(c^), 

up to e 2 , where the shear tensors are 

°ij = d(iVj) - ^-0, (cr^)^ = d(i<pdj)<p - ^00, (3.29) 

and = 5^diVj, 9^ = 5 13 di<f)dj<ft ■ Our aim is to obtain the forced NS equations, which 
turn out to be the constraint equations of the gravitational field equations ( |3.4|) at order 
e 3 . The solutions up to e 2 are enough to provide the forced NS equations, because higher 
order corrections do not make contribution in this order. With the solutions Q3.28 ), 
the constraint equations at order e 2 give the incompressible condition co diV l = 0, and 
at order e 3 give us with 

co {d T Vi + v j d jVi + d t P) - Vo d 2 v t = + fl q) . (3.30) 

These equations are just the temporal and spatial components of the equations in (|3.24|) 
up to e 3 . Here the external forces only have the spatial components as 

/<*> = - Xo (vidrfdtj)) + £ ($^84) , = -n (Ef 9) + tfB<j?>) , (3.31) 

where Ef 9 ^ = and Bf- 9 ^ = F^ are the background electric and magnetic fields 
respectively. Further we define the dynamic viscosity and the normalized forces as 

"M = W"o, f{% = fP/uo, /S = fi q) /^, (3.32) 
the forced incompressible NS equations then become 

d T Vi + v j d jVi + d t P - u [u] d\ = /gj + d lV l = 0. (3.33) 

If we consider the characteristic scale L ~ e _1 and the velocity v = V v i v% ~ the 
Reynolds number of the dual fluid 1Z e = vL/vla ~ A'kTjj [1 + (n fi )/ (T H so)}. 
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In addition, we can divide the electromagnetic forces into two parts // = + 
f- ex \ where 

ft x) = - n (Ef x) + xfiB$X) , ft d) = borwofrvi + v^d jVi + d t P), (3.34) 

where E^' = F^f 1 and = F^. Then is just the Lorentz force from the 

extra electromagnetic field, while from the boundary chemical potential. We now 
move the term to the left hand side of the NS equations. Defining the dynamic 
shear viscosity and the force density as 

v , , = = 1 f(f) = J±- f(f ) = li (3.35) 

we can rewrite the incompressible charged NS equations as 

d tVi + v j d jVi + d,P - u [s] d\ = f$ + fffi, dtf = 0. (3.36) 

In this case, the Reynolds number becomes lZ e = vL/u[ a ] oc T#, proportional to the 
temperature of the fluid. 



4. Holographic charged fluid at cutoff surface 

In this section, we will generalize the previous discussions to the case of dual fluid at a 
finite cutoff surface by using the method which is introduced in section.^. In this case, 
we only need to substitute 

9tt(r) = g-\r) = r 2 V(r)/f, g xx (r) = r 2 /f, (4.1) 

into the generic metric (|2.1|) , and introduce a finite cutoff surface in the Einstein- 
Maxwell system with a dilaton field and consider the non-relativistic expansions 

«9 r ~e°, dt ~ 0* ~ d i( j) ~ e\ <9 T ~P~e 2 . (4.2) 

The equations of motion to be solved are given by ( |3.2j ) and (|3.4|). 

4.1 Charged fluid in non-relativistic limit 

Since we are considering the bulk solution with a finite cutoff in this section, following 
31~], we choose the gauge that g rr = 0, g ra oc u a and g£ = 0. The perturbed solution 
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up to order e 2 with a Dirichlet boundary condition at the cutoff surface turns out to be 



d~s 2 = - 



+ 



1 + H(r)(r% 



,~2 



u a (x)dx a dr + — 



V(r) _ 



V(r 



—U a (x)u b (x) + Vab(x) 



doc doc 



nr){r%)^ + 2r c T(r)y/V^j a ab + r c 2 ^(r)(^) a6 



V(r c 



dx a dx b 



A =A a {r, x)dx a = A^ d \r, x) + A^ x \x) 



dx a + 



/x Q(r) r£ 1 , 



(4.3) 



T {r%)u a {x)dx a , 



$ =0(F) + r c ^(r) (V4</>) + r c 2 J^(r) (>4< 



where the gauge field relates to (|3.f 2|) through 



t U{r) 



L=u a (x), A^\x) 



A { r\x)„ 



u a (x) + ^A { r\x)^Ax). 



(4.4) 



and n l a {x) are given in (|2.27|) with a coordinate dependent velocity. Again keep in mind 
that we have already taken the non-relativistic limit ( |4.2| ) of the solution ( |4.3| ) up to 
order e 2 and higher order terms do not make contribution in the following calculations, 
they are reserved to only keep a covariant form. For the shear perturbations due to the 
boost, 



^y^dy, o ah = d {a u b) - —9. 



(4.5) 



and 9 = rj ab d a u b = V ab d a u b . For the shear perturbations due to the dilaton field 



p — 1 



2 /-r c / p-l _ p-1n _ -p 

7 /; (^)o6 = 0«, (4-6) 



y p+2 ^(y) 



where 6^ = V ab (d a <f)d b <j)) . The corresponding scalar perturbation equations give the 
following solutions 



£(r) 



/i(r c 



4p(p - l)r 2 0^7) 



Q(r) 



a r. 



4p(p - l)r 2 V r, 



r p-i 



-P-i / ' 



2p(p — l)r 2 



(4.7) 



Here two of the integration constants in Q4.7 ) have been determined via the Dirichlet 
boundary condition of /C(r) and Q(r). The other two integration constants a(r c ), h(r c ), 
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and the notation 

, N h(r c ) - a(r c ) r 2 q 2 £ 2 ( r p+1 \ 

g{r) = A 7mr^ (48) 

will be determined via choosing the Landau gauge of the stress energy tensor of dual 
fluid and the induced charge current, whose expressions are given by 

T ab (r c ) = 2(K Vab - K ab ) + lf\ Jj?\r c ) = N c F a c . (4.9) 

The counter term at the cutoff surface is 

Tj b ct) = -c(r c ) fe afe + ( 4040 - %(<90) 2 H , (4.10) 



where c(r c ) could be an arbitary function of r c . In fact, at the finite cuoff surface, the 
counter term is not necessary. To match the results at the AdS boundary in the case 
without the cutoff surafce, we here add the counter term and can take c(r c ) = 1. 

The stress energy tensor and current in the derivative expansion of the fluid dual 
to the perturbed solution ( (4.3| ) can be written as 

fabij c ) = t} b ] + + t} b 2) + ..., r} b 0) = u(r c )u a u b + p(r c )^, 
#(r c ) = + + + = n(r c )u a , 1 ' 1 

where 0(r c ) = p(r c ) +p(r c ). The zeroth order parts in derivative expansion are just 
the stress energy tensor and current of ideal charged fluid. The energy density p(r c ), 
pressure p(r c ), and charge density n(r c ) up to order e 2 in the non-relativistic expansion 
are given by 



P( r c) = P {r c ) ~ {r c P)p {r c ), p (r c ) = + — 



" 1 



r c V'(r c ) 

p{r c ) = p (r c ) - (r c P)p' o (r c ), p (r c ) = A== - Po{r c ), (4.12) 

W V i r c) 

(r c ) = n (r c ) - (r c P)n' Q (r c ), n (r c ) = U'(r c ). 



n 



The higher order terms in derivative expansion ( [4.1 1[) are the dissipative parts of 
the stress energy tensor and current, we denote them by T^ b %ssS> and Jat %SS ' '• Then using 
the Landau gauge up to e 2 , 

u*Tf ss) = =► h(r e ) = c(r e ), u a J^ ss) = =► a(r c ) = h(r c ), (4.13) 

so that q(r) = 0, and when r c — > oo we can recover the results ( gig) at the AdS 
boundary via transforming them to the coordinates in the boundary. 
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After imposing the Landau gauge (|4.13|) up to order e 2 , we have Ja = Ja = 



and 

'W = - 2 V ( r c)°ab, t} b 2) = -2r^(r c )(c^) ab - <^(r c )^P ab , 
where the cutoff dependent shear viscosity and entropy density are 



1 < 



4G p+2 rl 

And the coefficients associated with dilaton cf) are 



5 (r c ) 4tt' 



(4.14) 



(4.15) 



i ~t \ /t77 — V c(r c )r c F'(^c) 



(4.16) 



On the other hand, at the zeroth order, from ( |2.7|) and ( |4.12| ), one has 



^o(r c ) 



4vr£r c ^(^)' 



(4.17) 



as well as ^ (r c ) in fl4.12|) , and s (r c ) in (|4.15|) , thus one can show that the following 
thermodynamical relation still holds for the fluid at the cutoff surface 



T {r c )s (r c ) = Q {r c ) - h {r c )fl {r c ), £ (r c ) 



e[U(r c )-U(r h )] 
r c \/V{r c ) 



(4.18) 



where the chemical potential fi (r c ) is defined as the difference of the gauge potential 
between the horizon and the cutoff surface. The dimensionless coordinate invariant 
diffusivity defined in [[L4| is 



D c {r c ) = f {r c )^ Tr 



1 + ^o( r c)/i ( r c 



-1 



(4.19) 



£ (r c ) 4vr [ T (r c )s (r c ) 
which shows the dependence of the cutoff surface when the chemical potential is present. 

4.2 Forced Naiver-Stokes equations 



As what has been done in section |3.1| , we can redefine the stress energy tensor of the 
dual fluid as 



= u(r c )u a u b +p{r c )r] ab - 2fj (r c )a, 



ab j 



(4.20) 
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by moving out the part of dilaton field. By using the constraint equations at the cutoff 
surface for gravitational field ( [3.2 ip and Maxwell field, up to e 3 we can get the forced 
NS equations and the conservation equation of charge current 



?(<?) 



(4.21) 



where the external force from the dilaton field is 



ft"\r c ) = -Ur c ) [u a d a 4dM + £ (r c ) ( &>d a 



X { r c) = Vo{ r c) = ~p, 
in. 



i ( r c 



(p-2)£ 



Pi? ~ l )\/ V ( r c) 



p-1 



+ 



(p-2) 



c(r c )r c V'(r c ) 



i4.22) 



And the Lorentz force of the charged fluid is 



b ' 



, (ex) _ _ ,ja rp(ex) 
Jb ~ >- J (q) r ab ' 



(bd) 



~>- J (q) r ab ' 



(4.23) 



with the external Maxwell tensor = 2d[ a A^ x \r c ), and the induced Maxwell tensor 
Pf£f* = 2d[ a A ( j\ d \r c ), where the gauge potential on the cutoff surface is 



£U(r c 



(4.24) 



If we furthermore redefine the following stress energy tensor by moving out the terms 
associated with the boundary chemical potential 



Tab = T {r c )s {r c )u a u b + \p (r c ) - n Q {r c )fx {r c )} fj ab - 2fj (r c )a ab . 



(4.25) 



where f c = r c (l — P), when b = 1, one has U(r h ) = 0, and ji (b) {r c ) is just the chemical 
potential fl {r c ) of the fluid at the cutoff surface, the following relation holds up to 
order e 3 , 



(bd) bp=l oa-T-(s) _ f(4>) , f{ex) 
f ° 1 ab — Jb "T Jb ■ 



(4.26) 



In this case, the external forces only come from the dilaton field and the external 
electromagnetic field. To see the above results more clearly, we write the them up to 
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the desired order. For example, the solutions 



ds 2 = + 



2(: 



2£ 



-drdr + 



+ 



Yrs. dr 2 + 5 ij dx i dx j 
V (r c ) 



up to the second order e 2 are, 
r 2 f. V (r 



1 



2(3idx l dT 



r c \ \ r, 

,2 



-fadx^r + \P + 



r c V'{r c ) - (3 2 
2V(r c ) + 2 



+ 



r 2 V(r) ( rV'if) r c V'(r c ) 



r 2 c V(r c ) V V(r) V(r c ) 



Pdr 2 + 



1 - 



V(r c 

drdr + T-L(r)(r 2 9^)dTdr 
V(r 



V(r c ] 



{fdr 2 + PiPjdx'dx*) 



+ 



A 



r 2 /C(r) 
r 2 F(r c ) 
£U{r) 



r%)drdr + — 2r c F(r)^V(r c )a lJ + r 2 F^r)(a 4 



2 



t/(r) 2V(r c ) 
HoQ{r) r p ~\ 



P 



dr 



dx { dx j , ( 4 - 27 ) 
£U(r) 



.1 V C 



$ = 

and the shear components are given by 

~ 5 
= du(3j) — —d k (3 k , 



V 



6 Hh 

— 6 4 
P 



(4.28) 



where 6 = d^&iPj, 9$ = 5^di4>dj(p. At order e 2 , the temporal component of the forced 
NS equations in ( 4.21|) turns out to be 



w (r c )^ = 0, tu (r c ) = r c V'{r c )/{l^/V{^)), (4.29) 
which leads to the incompressible condition, and the spatial component at order e 3 is 
w (r c ) (5 T 0i + /S^ft) - r^ir^P - Vo (r c )d 2 (3 t = fi*\r c ) + fi q \r c ), (4.30) 
where only the spatial components of external forces remain 
jf } (r c ) = " X (r c ) (Fdrfdrf) + | (r c ) (& fdfd^] 



ft\r B )= -n (r c ) ffi+pF™ 



(bd) 



If further redefine 



^o( r c) 



^o( r c) 



(4.31) 
(4.32) 

(4.33) 
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we can obtain the incompressible forced NS equations as 



d T fr + Pdtfi + diP M - D H d 2 fr = + /g, c^ 1 = 0. (4.34) 



Consider the characteristic scale of perturbations L ~ e 1 and velocity /3 = \f^fi % ~ e, 
the Reynolds number of the fluid turns out to be 



f r ) = JL. K _J_ = = 4vrf fr ) 



x n {r c )fj, (r c ] 
T (r c )s (r c ) 



(4.35) 



We can see that for uncharged black brane where the chemical potential vanishes, 
the Reynolds number of the dual fluid is proportional to the local temperature at 
the cutoff surface. Thus, when the cutoff surface approaches the horizon, the local 
temperature as well as the Reynolds number become larger and larger, and the fluid 
may become unstable. This instability may relate to the superluminal hydrodynamic 
sound modes when the cutoff surface is sufficiently close to the horizon []3T| , f49| . In 
addition, as in section |3.2| , we can divide the electromagnetic forces into two parts: 
R g) (r c ) = fT\r c ) + fT\r c ), where 

ft\r c ) = -n (r c ) (d T Af x) - ~d,A^ + vW 3 AJ x) - JdJ^) , (4.36) 
ft d \r c ) = n (r c ) [fi (r c ) (d T f3 { + - r^r^P] . (4.37) 

Define 

5 _ n (r c )r c fx' (r c ) - r c p' (r c ) ~ m _ f^\r c ) ~ {ex) _ f^\r c ) 



f (r c )s (r c ) ' [sl1 f (r c )s (r c y [sl1 T (r c )s (r c ) ' 

(4.38) 

and the dynamical shear viscosity z>[ s ] = r) (r c ) /[T (r c )s (r c )], the forced NS equations 
(|4.34j) then become 



d T 0i + Fdrfi + dii\ a] - u [s] d 2 fr = ft® + f™, = 0. (4.39) 

In this case, the external forces only come from the additional electromagnetic field and 
dilaton field. And the Reynolds number has the form TZ e (r c ) oc v^{r c ) = 47rT (r c ). 

5. Conclusions 

We have studied the thermodynamics and non-relativistic hydrodynamics of the holo- 
graphic fluid at a finite cutoff surface. As a calculation example, we have considered the 
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case with the Gauss-Bonnet gravity. The isentropic flow and shear viscosity of the dual 
fluid have been obtained. The radial Einstein equation implies the isentropy of RG flow 
was first proposed in [ 14] , and here we have generalized the discussion to the case of the 
Gauss-Bonnet gravity. The isentropy of RG flow can also be considered as an adiabatic 
process of the dual fluid. Note that instead of the entropy associated with the holog- 
raphy screen pD| |5T , here the total entropy of the dual fluid is the horizon entropy 
of the background black brane, hence cutoff independent. We have given a general 
formula (|2.50|) on the ratio of shear viscosity over entropy density of the fluid dual to 
the Gauss-Bonnet gravity. It shows that the ratio is independent of the cutoff surface. 
Namely it does not run with the cutoff surface. This may explain why the membrane 

1& 



paradigm |j2| |3[ and the AdS/CFT correspondence 
give the same result. 
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The main goal of this paper is to give the expressions of external force, coming from 
the bulk matters, of holographic fluid in the non-relativistic limit. To this end, we have 
considered an Einstein-Maxwell-dilaton system with a negative cosmo logical constant. 
By using the non-relativistic fluid expansion method, we have solved the system up 
to the second order of non-relativistic fluid expansion parameter e, and obtained the 
incompressible forced NS equations of the dual fluid at the AdS boundary and at a 
finite cutoff surface, respectively. The concrete expressions of external forces from the 
dilaton field and Maxwell field have been given. Here we have taken an new scaling of 
the dilaton field di<p ~ e 1 in the non-relativistic limit, so that the external force provided 
by the dilaton field fl ~ e 3 , meeting the scaling symmetry of the forced NS equations. 
Actually, in the derivative expansion of the stress energy tensor, these terms such as 
di(pdj(f> are the second order dissipative terms of the dual fluid [fL2"l , and such terms 
may appear in the superfluid components f62| . However, in the non-relativistic limit, 
we move them to the right hand side of the NS equations as external force terms. Note 
that here we have considered the case of minimal coupling of dilaton field. It would be 



interesting to extend this study to the case with non- minimal coupling [53, 54 



It turns out the Reynolds number of dual fluid is proportional to the local tem- 
perature of the cutoff surface in the uncharged case. Thus when the cutoff surface 
approaches to the event horizon of the black brane background, the local temperature 
and thus the Reynolds number become larger and larger. This indicates that the dual 
fluid will become unstable when the cutoff surface is close enough to the event horizon. 
It would be interesting to further study the stationary turbulence by using the forced 
NS equations derived in this paper. 
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